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ABSTRACT 

We identify SDSS 011009.09+132616.1, SDSS 030308.35+005444.1, 

SDSS 143547.87+373338.5 and SDSS 154846.00+405728.8 as four eclipsing white 
dwarf plus main sequence (WDMS) binaries from the Sloan Digital Sky Survey, and 
report on follow-up observations of these systems. SDSS 0110+1326, SDSS 1435+3733 
and SDSS 1548+4057 contain DA white dwarfs, while SDSS 0303+0054 contains a cool 
DC white dwarf. Orbital periods and ephemerides have been established from multi- 
season photometry. SDSS 1435+3733, with Forb = 3h has the shortest orbital period 
of all known eclipsing WDMS binaries. As for the other systems, SDSS 0110+1326 
has Porb = 8h, SDSS 0303+0054 has Porb = 3.2 h and SDSS 1548+4057 has 
-forb = 4.4 h. Time-resolved spectroscopic observations have been obtained and 
the Ha and Can AA 8498.02,8542.09,8662.14 triplet emission lines, as well as the 
NaiAA 8183.27,8194.81 absorption doublet were used to measure the radial velocities 
of the secondary stars in all four systems. A spectral decomposition/fitting technique 
was then employed to isolate the contribution of each of the components to the 
total spectrum, and to determine the white dwarf effective temperatures and surface 
gravities, as well as the spectral types of the companion stars. We used a light curve 
modelling code for close binary systems to fit the eclipse profiles and the ellipsoidal 
modulation/reflection effect in the light curves, to further constrain the masses and 
radii of the components in all systems. All three DA white dwarfs have masses of 
-^WD ~ 0.4 — 0.6 M0, in line with the expectations from close binary evolution. 
The DC white dwarf in SDSS 0303+0054 has a mass of Mwd > 0.85 Mq, making it 
unusually massive for a post-common envelope system. The companion stars in all 
four systems are M-dwarfs of spectral type M4 and later. Our new additions raise 
the number of known eclipsing WDMS binaries to fourteen, and we find that the 
average white dwarf mass in this sample is < A/wd >= 0.57 ± 0.16 Mq, only sHghtly 
lower than the average mass of single white dwarfs. The majority of all eclipsing 
WDMS binaries contain low-mass (< 0.6 M©) secondary stars, and will eventually 
provide valuable observational input for the calibration of the mass-radius relations 
of low-mass main sequence stars and of white dwarfs. 

Key words: binaries: close - binaries: eclipsing - stars: fundamental 
parameters - stars: late-type - stars: individual: SDSS 011009.09+132616.1, 
SDSS 030308.35+005444.1, SDSS 143547.87+373338.5, SDSS 154846.00+405728.8 - 
white dwarfs 
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1 INTRODUCTION 

White dwarfs and low-mass stars represent the most com- 
mon types of steUar remnant and main sequence star, re- 
spectively, encountered in our Galaxy. Yet, despite be- 
ing very common, very few white dwarfs and low mass 
stars have accurately determined radii and masses. Con- 
sequently, the fini te temperat u re mass-radius re lation of 
white dwarfs (e.g. I WoodI fl995l : IPanei et al.l 200C ) rern ains 
largely untested by observations ( Provencal et al. Il998l ). In 
the case of low mass stars, the empirical measurements con- 
sistently result in radii up to 15% larger and effective tem- 
peratures 400K or more below the values p redicted by the- 
ory (e.g. iRibasI l2006l : lL6pez-Mora"i5 I2007D . This is most 
clearly demon strated using low-mass eclipsing binary stars 
Bavless fc Orosz .2006'). but is also pr esent in field stars 



Table 1. Full SDSS names and u,g,r,i,z magnitudes 
of SDSS 0110+1326, SDSS 0303+0054, SDSS 1435+3733 and 
SDSS 1548+4057. 



Berger et al]|2006l : [Morales et"aLll2008l) and the host stars 



of transiting extra-solar planets (Torres l2007f ). 

Eclipsing binaries are the key to determine ac- 
curate steUar masses and radii (e.g. lAndersenI Il99ll : 
ISouthworth fc ClausenlboOTp . However, because of their in- 
trinsic faintness, very few binaries containing white dwarfs 
and/or low mass stars are currently known. Here, we report 
the first results of a programme aimed at the identification of 
eclipsing white dwarf plus main sequence (WDMS) binaries, 
which will provide accurate empirical masses and radii for 
both types of stars. Because short-period WDMS binaries 
underwent common envelope evolution, they are expected 
to contain a wide range of white dwarf masses, which will 
be important for populating the empirical white dwarf mass- 
radius relation. Eclipsing WDMS binaries will also be of key 
importance in filling in the mass-radius relation of low-mass 
stars at masses < 0.6 Mq. 

The structure of the paper is as follows. The target se- 
lection for this programme is described in Sect. (2] In Sect.[3] 
we present our observations and data reduction in detail. We 
determine the orbital periods and ephemerides of the four 
eclipsing WDMS binaries in Sect.]!] and measure the radial 
velocities of the secondary stars in Sect.O In Sect. |B] we de- 
rive initial estimates of the stellar parameters from fitting 
the SDSS spectroscopy of our targets. Basic equations for 
the following analysis are outlined in Sect.[71 In Sect.[S]we 
describe our fits to the observed light curves, and present our 
results in Sect.[9l The past and future evolution of the four 
stars is explored in Sect. 1101 Finally, we discuss and sum- 
marise our findings, including an outlook on future work in 
Sect.[lTl 



2 TARGET SELECTION 

We have selected eclipsing SDSS WDMS binaries based on 
the available information on the radial velocities of their 
companion stars, and/or evidence of a strong reflection ef- 
fect. 

Initially, we used SDSS spectroscopy to mea- 
sure the radial velocity of the companion star either 
from the Nai AA 8183.27,8194.81 absorption doublet, or 
from the Ha emission line (see iRebassa-Mansergas et ahl 
l2007l for details). SDSS J030308.35+005444.1 (henceforth 
SDSS 0303+0054) and SDSS J143547.87+373338.6 (hence- 
forth SDSS 1435+3733) exhibited the largest secondary star 
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radial velocities among ~ 1150 WDMS binaries which have 
SDSS spectra of sufficiently good quality, 287 km s^^ and 
335kms"\ respectively. For SDSS J011009. 09+132616.1 
(henceforth SDSS 0110+1326) two SDSS spectra are avail- 
able, which differ substantially in the strength of the emis- 
sion lines from the heated companion star. Photometric time 
series (Sect.[3ll revealed white dwarf eclipses in all three ob- 
jects. SDSS 1435+3733 has bee n independently identified as 
an eclipsing WDMS binary bv lSteinfadt et all (120081 ). 

As the number of SDSS WDMS binaries with 
known orbital periods, Porb, and radial velocity ampli- 
tudes, Ksec, is steadily growi ng (ISchreiber et al.l l2008l : 
iRebassa-Mansergas et ahl l2008bl ). we are now in the po- 
sition to further pinpoint the selection of candidates for 
eclipses: with Porb and Ksec from our time-series spec- 
troscopy, and MwD and Msec fro m our spectral decomposi- 
tion/ fitting of the SDSS spectra (|Rebassa-Mansergas et aU 
[2003), we can estimate the binary inclination from Kepler's 
third law. In the case of SDSS J154846.00+405728.8 (hence- 
forth SDSS 1548+4057), the available information suggested 
i 85°, and time-series photometry confirmed the high in- 
clination through the detection of eclipses. 

Full coordinates and SDSS u, g, r, i, z point-spread func- 
tion magnitudes of the four systems are given in Table [1] 



3 OBSERVATIONS AND DATA REDUCTION 

Follow-up observations of the four systems - time series CCD 
photometry and time-resolved spectroscopy - were obtained 
at five difi'erent telescopes, namely the 4.2m William Her- 
schel Telescope (WHT), the 2.5m Nordic Optical Telescope 
(NOT), the 2.2m telescope in Calar Aho (CA2.2), the 1.2m 
Mercator telescope (MER) and the lAC 0.8m (IAC80). A 
log of the observations is given in Table[51 while Fig.[T]shows 
phase-folded light curves and radial velocity curves of all the 
systems. A brief account of the used instrumentation and the 
data reduction procedures is given below. 



3.1 Photometry 

Photometric observations of the four targets were obtained 
at all five telescopes. In every telescope set-up, care has been 
taken to ensure that at least 3 good comparison stars were 
available in the science images, especially in the cases when 
the CCDs were windowed. 



3.1.1 WHT 4.2m 

The observations were carried out using the AUX-port im- 
ager, equipped with the default 2148 x 4200 pixel E2V 



Four eclipsing white dwarf main sequence binaries 3 



Table 2. Log of the observations. 
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CCD44-82 detector, with an unvignetted, circular field-of- 
view (FOV) of 2.2' in diameter. For the observations, the 
CCD was binned (4 x 4) to reduce readout time to 4 s. 
SDSS 1548+4057 was observed with a Johnson R fiher. The 
images were de-biased and flat-fielded within MfDAS and 
aperture photometry was carried out using SEXTRACTOR 
l|Bertin fc Arnout"3 1 19961 ) . A full account of the em ployed 
reduction pipeline is given bv lCansicke et all (|2004l ). 
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Figure 1. Phase-folded light- and radial velocity curves of the 
four systems. From top to bottom (two panels for each sys- 
tem): CA2.2 /-band light curve and Call radial velocity curve 
of SDSS 0110+1326, CA2.2 fitlerless hght curve and Nai radial 
velocity curve of SDSS 0303+0054, IAC80 7-band light curve and 
Nal radial velocity curve of SDSS 1435+3733 and WHT ij-band 
light curve and Nal radial velocity curve of SDSS 1548+4057. 



3.1.2 NOT 2.5m 

The observations were carried out using the Andalucia Faint 
Object Spectrograph and Camera (ALFOSC), equipped 
with a 2k X 2k pixel E2V CCD42-40 chip, with a FOV 
of 6.5' X 6.5'. The CCD was binned (2 x 2) to reduce 
readout time. Filterless observations were obtained for 
SDSS 1548+4057. The reduction procedure for the NOT 
data was the same as the one used for the WHT data. 



3.1.3 Calar Alto 2.2m 

The observations were carried out using the Calar Alto 
Faint Object Spectrograph (CAFOS) equipped with the 
standard 2k x 2k pixel SITe CCD (FOV 16' x 16'). For 
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the observations, the CCD was windowed to reduce readout 
time to 10 s. Observations were carried out using Johnson 
V, R, I, and "Roser" BvEl filters, as well as without any fil- 
ter. In detail, we used / and BV for SDSS 0110+1326, R 
and BV for SDSS 0303+0054, and V, BV, and no filter for 
SDSS 1435+3733. The reduction procedure was the same as 
described above. 



3.1.4 Mercator 1.2m 

Filterless photometric observations have been obtained for 
SDSS 0110+1326 and SDSS 0303+0054 using the MERca- 
tor Optical Photometric imagEr (MEROPE) equipped with 
a 2k X 2k EEV CCD chip (FOV 6.5' x 6.5'). For the obser- 
vations, the CCD was binned (3 x 3) to reduce readout time 
to about 8 s. The reduction procedure for the Mercator data 
was the same as above. 



3.1.5 lACO.Sm 

Johnson J-band photometry has been obtained at 
the lACO.Sm telescope for SDSS 0110+1326 and 
SDSS 1435+3733, while V- and i?-band photometry 
has been obtained for SDSS 1548+4057. The telescope 
was equipped with the SI 2k x 2k E2V CCD (FOV 
10.25' x 10.25'), which was binned (2 x 2) and windowed to 
reduce readout time to lis. Data reduction was performed 
within IRAF0. After bias and flat-field corrections the 
images were aligned and instrumental magnitudes were 
measured by means of aperture photometry. 



3.2 Spectroscopy 

All spectroscopy presented in this work was acquired using 
the William Herschel Telescope (WHT) and the Interme- 
diate dispersion Spectrograph and Imaging System (ISIS). 
For the red arm the grating used was either R316R or 
R600R, and for the blue arm the grating was either R600B or 
R1200B. For all observations the blue-arm detector was an 
EEV 2kx4k CCD. In 2006 July and September the red-arm 
detector was a Marconi 2kx4k CCD, and in later observ- 
ing runs a high-efficiency RED+ 2kx4k CCD was used. In 
all cases the CCDs were binned spectrally by a factor of 2 
and spatially by factors of 2-4, to reduce readout noise, and 
windowed in the spatial direction to decrease the readout 
time. 

Data reductio n was undertaken using optimal extrac- 
tion (iHornel Il98^ ) as implemented in the pamel43 code 
iMarsh 1983), which also makes use of the STARLINPO pack- 
ages FIGARO and KAPPA. Telluric lines removal and flux- 



1 A BG39/3 filter centred on 4977 A with a full width at half 
maximum (FWHM) of 1559 A 

^ IRAF is distributed by the National Optical Astronomy Obser- 
vatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under contract with the National 
Science Foundation, |http://iraf.noao.edu l 

3 PAMELA and MOLLY were written by TRM and can be found at 
|http : //www. Warwick, ac .uk/go/trmeirsh 
* The Starlink Software can 
|http : //starllnk. jach.hawali ■ edu/| 



calibration was performed separately for each night, using 
observations of flux standard stars. 

Wavelength calibrations were obtained in a standard 
fashion using spectra of copper-argon and copper-neon arc 
lamps. For SDSS 0303, arc lamp exposures were obtained 
during the spectroscopic observations and the wavelength 
solutions were interpolated from the two arc spectra brack- 
eting each spectrum. For the other objects we did not obtain 
dedicated arc spectra, to increase the time efficiency of the 
observations. The spectra were wavelength-calibrated using 
arc exposures taken at the beginning of each night, and drift 
in the wavelength solution was removed using measurements 
of the positions of the A7913 and A6300 night-sky emission 
lines. This procedure has been found to work well for sp ec- 
tra at red wavelengths (see lSouthworth et al.ll200^ . l2008l ). In 
the blue arm, a reliable correction of the wavelength zero- 
point was not possible as only one (moreover weak) sky-line 
(Hgl A4358) is available, and consequently, these spectra 
are not suitable for velocity measurements. The reciprocal 
dispersion and resolution for the R316R grating are approx- 
imately 1.7 Apx~^ and 3.3 A, and for the R600R grating are 
0.89 Apx~^ and 1.5 A, respectively. 



4 ORBITAL PERIODS AND EPHEMERIDES 

Mid-eclipse times for all four systems were measured from 
their light curves as follows. The observed eclipse profile was 
mirrored in time around an estimate of the eclipse centre. 
The mirrored profile was then overplotted on the original 
eclipse profile and shifted against it, until the best over- 
lap between the points during ingress and egress was found. 
Given the sharpness of eclipse in/egress features, we (con- 
servatively) estimate the error in the mid-eclipse times to be 
comparable to the duty cycle (exposure plus readout time) 
of the observations. Table [3] lists the mid-eclipse times. An 
initial estimate of the cycle count was then obtained by fit- 
ting eclipse phases (0g''=<="'™d _ ^&t-^~2 ^^^^ ,^ wide range 
of trial periods. Once an unambiguous cycle count was es- 
tablished, a linear eclipse ephemeris was fitted to the times 
of mid-eclipse. For SDSS 1435 +3733 we have a l so us ed the 
mid-eclipse times provided by iSteinfadt et al.l l|2008l ) . The 
resulting ephemerides, with numbers in parenthesis indicat- 
ing the error on the last digit, are 



be 



found 



at 



ro(HJD) = 2453994.44787(9) + 0.3326873(1)E (1) 

for SDSS 0110+1326, that is, Porb = 7.984495(3)h 

ro(HJD) = 2453991.6164(1) + 0.13443772(7)E (2) 

for SDSS 0303+0054, that is, Porb = 3.226505(l)h 

To(HJD) = 2454148.70361(6) + 0.1256311(1)E (3) 

for SDSS 1435+3733, that is, Poa = 3.015144(2)h and 

ro(HJD) = 2454592.57135(6) + 0.18551774(4)E (4) 

for SDSS 1548+4057, that is, Porb = 4.45242576(4)h. 

These ephemerides were then used to fold both the pho- 
tometric and the spectroscopic data over phase (Fig.[l}. 
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Table 3. Mid-eclipse timings, cycle number and the difference between observed and computed eclipse times using the 
ephemerides provided in equations IH4I The large O — C values for the first five eclipses of SDSS 1548+4057 are due to the poor 
time resolution of the lAC 80 light curves (Table[2ll. Mid-eclipse times for SDSS 1435+3733 from lSteinfadt et al.l (|2008|'| are also 
included. 
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Figure 2. SDSS spectra of the four systems. Top to bottom: SDSS 0110+1326, SDSS 1435+3733 SDSS 1548+4057 and SDSS 0303+0054. 
The radial velocity of the companion star in SDSSOllO+1326 was measured from the Ha and Call emission lines, whereas for the other 
three systems, we used the Nal absorption doublet 



5 RADIAL VELOCITIES OF THE 
SECONDARY STARS 

Figure [2] shows the SDSS spectra of our four targets. Mea- 
suring radial velocity variations of the stellar components 
needs sharp spectral features. Determining the radial ve- 
locity amplitude of the white dwarf, Kwd, is notoriously 
difficult because of the width of the Balmer lines, and our 
spectroscopic data are insufficient in quantity and signal-to- 



noise ratio for the approach outlined e.g. bv lMaxted et al] 
(|2004i fl. We hence restrict our analysis to the measurement 
of the radial velocity amplitude of the secondary star, i^scc- 



^ Ultraviolet spectroscopy allows accurate measurements of 
^WD using narro w metal li nes originatin g in the white 
dwarf photosphere , l lO'Brien et a l. 2001; O'Don oghue et al.ll2003l ; 
iKawka et al.ll2007l) . but at the cost of space- based observations. 
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Figure 3. Sample line profile fits (gray) to the observed spectra 
(black) of SDSS 0303+0054 (left panel, Na I AA 8183.27,8194.81 
absorption doublet) and of SDSS 0110+1326 (right panel. Call 
emission triplet). The rest wavelengths are indicated by the dot- 
ted vertical lines. The velocity scale at the top is relative to the 
blue-most component of each multiplet. The top spectra are from 
SDSS, the middle and bottom spectra were obtained at the WHT. 
Radial velocity variations in these systems are obvious to the eye. 

5.1 SDSS 0110+1326 

The spectrum of SDSS 0110+1326 displays the 
Ca 11 AA 8498.02, 8542.09,8662.14 emission triplet as 
well as an Ha emission line. We determined the velocities of 
the Ha emission line by fitting a second-order polynomial 
plus a Gaussian emission line to the spectra. For the Can 
triplet, we fitted a second-order polynomial plus three 
Gaussian emission lines with identical width and whose 
separations were fixed to the corresponding laboratory 
values (see Fig.|31 right panel). The Ha and Can radial 
velocities were then separately phase-folded using the 
ephemeris Eq. and fitted with a sine wave. The phasing 
of the radial velocity curves agreed with that expected 
from an eclipsing binary (i.e. red-to-blue crossing at orbital 
phase zero) within the errors. The resulting radial velocity 
amplitudes are Knee ,Ha = 200.1±4.8 kms~^ with a systemic 
velocity of 'yua = 19.4 ± 4.1kms~^ for the Ha line; and 
J^aec.Caii = 178.8 ± 2.4 km s" ^ and7caii = 15.2± 2.4 kms"^ 
for the Can line (see also Fig.[T]). 

We decided, following a suggestion by the referee,to in- 
vestigate whether the emission lines originate predominantly 
on the illuminated hemisphere of the secondary star. If that 
is the case - keeping in mind that we see the system al- 
most edge-on - we expect a significant variation of the line 
strength with orbital phase, reaching a maximum around 
phase (j> — 0.5 (superior conjunction of the secondary) and 
almost disappearing around (fi — 0.0. Figure|4] shows average 
spectra of SDSS 0110+1326, focused on the Ha line and the 
Can triplet for — 0.0 and (f} = 0.5. It is apparent that the 
Ha and Can emission lines are very strong near (j) = 0.5. The 
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Figure 4. Average spectra of SDSS 0110+1326. Left panel: Ha 
line, right panel: Can triplet. The binary phase of each spectrum 
is clearly marked. Prior to averaging, each spectrum was shifted 
to the restframe of the secondary, using the measured i^sec and 
7 values. For each line and binary phase, three spectra were av- 
eraged. 

Ha emission line is very weak near — 0.0, and Can is seen 
in absorption. The equivalent width (EW) of the blue-most 
component of the Ca n triplet is shown in Fig.[S] (top panel) 
as a function of orbital phase. Given the spectral resolution 
and quality of our data, measuring the EW of the Ha emis- 
sion line is prone to substantial uncertainties, as it is embed- 
ded in the broad Ha absorption of the white dwarf photo- 
sphere (see Fig. [2] again) , but it generally follows a simlilar 
pattern as the one seen in the Can triplet. This analysis 
supports our assumption that the emission lines originate 
on the irradiated, inner hemisphere of the secondary star. 
Hence, the centre of fight of the secondary star is displaced 
towards the Lagrangian point Li, with respect to the centre 
of mass. The emission lines trace, therefore, the movement 
of the centre of light and not the centre of mass and as a 
result, the Kscc values measured from either the Ha or the 
Can emission lines are very unlikely to represent the true 
radial velocity amplitude of the secondary star, but instead 
give a lower limit to it. 

The fact that ifsec measured from the Ha line is larger 
than that from the Can line suggests that the Ha emission 
is distributed in a slightly more homogeneous fashion on 
the secondary, and illustrates that correcting the observed 
Kacc for the effect of irradiation is not a trivial matter. We 
describe our approach to this issue in Sect. 18.21 

5.2 SDSS 0303+0054, SDSS 1435+3733, and 
SDSS 1548+4057 

The NaiAA 8183.27,8194.81 absorption doublet is a strong 
feature in the spectra of SDSS 0303+0054, SDSS 1435+3733, 
and SDSS 1548+4057. We measured the radial velocity vari- 
ation of the secondary star in SDSS 0303+0054 by fitting 
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Figure 5. Equivalent widths of the lines used to measure radial 
velocities in all four systems with respect to the orbital phase. 
Top panel: blue-most component of the Ca II emission triplet in 
SDSS 0110+1326. Lower three panels: blue-most component of 
the Na I absorption doublet in SDSS 0303+0054, SDSS 1435+3733 
and SDSS 1548+4057 respectively. The dashed lines indicate the 
mean EW. The strong variation of the EW of the Call emission 
line is obvious. No significant variation as a function of phase is 
observed for the other three systems. A full orbital cycle has been 
duplicated for clarity. 



Table 4. Summary of the radial velocity measurements for all 
four systems. 



System 


Line 


Kscc [kms ^] 


7 [kms ^] 


SDSS 0110+1326 


Ha 


200.1 + 4.8 


19.4 + 4.1 




Call 


178.8 + 2.4 


15.2 ± 2.4 


SDSS 0303+0054 


Nal 


339.7+ 1.9 


-4.0+ 1.4 


SDSS 1435+3733 


Nal 


260.9 + 2.9 


47.4 + 2.2 


SDSS 1548+4057 


Nal 


274.7+ 2.6 


-7.4+2.2 



this doublet with a second-order polynomial plus two Gaus- 
sian emission lines of common width and a separation 
fixed to the corresponding laboratory value (see Fig.O left 
panel). The same method was applied to the spectra of 
SDSS 1435-^3733 and SDSS 1548-^4057. A sine-fit to each 
of the radial velocities data sets, phase-folded using the 
ephemerides Eq. (Hjl-Q, gives the radial velocity of the sec- 
ondary star, Kscc, and the systemic velocity 7 for each sys- 
tem. The results of the radial velocity measurements are 
summarised in Table |4l 

No significant variation in the strength of the Nai dou- 
blet was observed as a function of orbital phase, and we 
hence assume that T^soc measured from this doublet reflects 
the true radial velocity amplitudes of the secondary stars in 
SDSS0303-H0054, SDSS 1435+3733, and SDSS 1548+4057. 
This is illustrated in Fig. [5] where we plot the equivalent 
widths against orbital phase for these three systems. 




4000 5000 6000 7000 BOOO 9000 

Wavelength [A] 

Figure 6. Two-component fit to the spectrum of 
SDSS 0110+1326. The top panel shows the spectrum of the 
object as a solid black line and the two templates, white dwarf 
and M-dwarf, as dotted lines. The bottom panel shows the 
residuals from the fit. 



6 SPECTROSCOPIC STELLAR PARAMETERS 

We used the spe ctral decomposition/fittin g met hod de- 
scribed in detail in iRebassa-Mansergas et al.l (|2007l ) to esti- 
mate the white dwarf effective temperatures (Tes) and sur- 
face gravities (logg), as well as the spectral types of the 
companion stars, for our four targets from their SDSS spec- 
troscopy. 

Briefiy, the method employed is the following: as a first 
step a two-component model is fitted to the WDMS binary 
spectrum using a grid of observed M-dwarf and white dwarf 
templates. This step determines the spectral type and fiux 
contribution of the M-dwarf component as shown in Fig.|6l 
After subtracting the best-fit M-dwarf, the residual white 
dwarf spectr um is fitted w i th a grid of white dwarf model 
spectra from lKoester et ahl (|2005l ). We fit the normalised H/3 
to He line profiles, omitting Ha which is most severely con- 
taminated by the continuum and/or Ha emission from the 
companion star. Balmer line profile fits can lead to degener- 
acy in the determination of Toft and log g, as their equivalent 
widths (EWs) go through a maximum at T^s — 13000 K, 
which means that fits of similar quality can be achieved for 
a "hot" and "cold" solution. In order to select the physically 
correct solution, we also fit the continuum plus Balmer lines 
over the range 3850-7150 A. The resulting T^a and logg are 
less accurate than those from line profile fits, but sensitive 
to the slope of the spectrum, and hence allow in most cases 
to break the degeneracy between the hot and cold solutions. 
Figure [7] illustrates this procedure. Once Tbs and logg are 
determined, the white dwarf mass and radius can be esti- 
mated using an updated version of iBergeron et all 's (1995) 
tables. The results, after applying this method to our four 
systems, are summarised in Table [S] The preferred solution 
("hot" or "cold") is highlighted in bold fond. 

In the case of SDSS 0303+0054, which contains a DC 
white dwarf (Fig. [2]), the spectral decomposition results in 
Sp(2) — M4.5 ± 0.5 for the secondary star. The subsequent 
fit to the residual white dwarf spectrum is not physical 
meaningful, as the white dwarf in SDSS 0303+0054 does 
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Figure 7. Spectral model fit to the white dwarf in 
SDSS 0110+1326, obtained after subtracting the best-fit M-dwarf 
template. Top left panel: best-fit (black lines) to the observed H/3 
to He (gray lines, top to bottom) line profiles. The model spectra 
and observations have been normalised in the same way. Top right 
panel: 1, 2 and 3 a contour plots in the T^ff-logg plane. The black 
contours refer to the best line profile fit, the red ones (which col- 
lapse into a dot on the scale of the plot) to the fit of the spectral 
range 3850-7150 A. The dashed line indicates the occurrence of 
maximum H/3 equivalent width. The best "hot" and "cold" line 
profile solutions are indicated by black dots, while the best fit 
to the whole spectrum by a red one. Bottom panel: the residual 
white dwarf spectrum resulting from the spectral decomposition 
and their fiux errors (gray lines) along with the best-fit white 
dwarf model (black line) in the 3850— 71 50 A wavelength range 
(top) and the residuals of the fit (gray line, bottom). 



not exhibit Balmer lines. The fit to the overall spectrum, 
which is purely sensitive to the continuum slope, suggests 
Toft < 8000 K, which is consistent with the DC classification 
of the white dwarf. 



7 BASIC EQUATIONS 

Here, we introduce the set of equations that we will use to 
constrain the stellar parameters of our four eclipsing WDMS 
binaries. In a binary system, where a white dwarf primary 
and a main-sequence companion, with masses Mwd and 
Msec respectively, orbit with a period Porb around their com- 
mon centre of mass at a separation a = awD + Oaec where 
awD Mwd ~ cisoc Msec, the orbital velocity K of either star, 
as observed at an inclination angle i, is 

Kj = ^^sini j = WD, sec (5) 

-forb 

assuming circular orbits. Given Kepler's law 

3 G (MwD + Msec) Port 



(6) 



Table 5. Summary of the results obtained for SDSS 0110-f 1326, 
SDSS 1435+3733 and SDSS 1548+4057 from our spectral decom- 
position technique. The preferred set of parameters for each sys- 
tem is highlighted in bold fond. 



Solution Mwd[M0] 


logs 


Tcft[K] 


Sp(2) 


SDSS 0110+1326 
Hot 0.47+0.02 

Cold 1.2 ± 0.03 


7.65+0.05 

9 + 0.04 


25891+427 

9619 + 23 


M4+1 

M4+ 1 


SDSS 1435+3733 
Hot 0.40 ± 0.05 
Cold 0.41+0.05 


7.58 + 0.11 
7.62+0.12 


12536 + 438 
12536+488 


M4.5 + 0.5 
M4.5+0.5 


SDSS 1548+4057 
Hot 0.43 + 0.16 
Cold 0.62+0.28 


7.64 + 0.31 
8.02+0.44 


14899 + 1300 
11699+820 


M6 + 0.5 
M6+0.5 



Mwd + Msec 
a = awD 

one can obtain the two mass functions 

I Kin 7^ 

/ (Msec) = 



/(Afv 



(Msec sini) 

(MwD + Msec 

(Mwd sin i) 



Porb A WD 



2tvG 



P uK'^ 

^ or b-* ^s( 

2tvG 



< Ms, 



< Mwd 



(7) 



(8) 



(9) 



(Mwd + Msec) 

which give strict lower limits on the masses of the compo- 
nents. From these two equations we get 

Msec KwD 



q = 



Mv 



ifs, 



(10) 



Therefore, the knowledge of the radial velocities of both stars 
can immediately yield the mass ratio q of the system, using 
Eq. (|10|l . In our case, since we lack a measurement for the 
radial velocity of the white dwarf Kwd, a more indirect 
approach needs to be followed. Re-arranging Eq. (|9]) for sini 
yields 



Porb K^c (Mwd + Ms cc 



-I 1/3 



2nG 



^WD 



whereas re-arranging Eq. (|9]) for Ks, 



2nG sin^ 



J" WD 



orb (A/wD + Afsec) 



■ yields 

1/3 



(11) 



(12) 



and using 



We make use of Equations (llip and (|12p in the light 
curve fitting process described in the next section. 

For the DC white dwarf in SDSS 0303+0054, where we 
lack a spectroscopic estimate of Mwd, we can use Eq. © 
to get a rough, first estimate of the white dwarf mass. In 
Fig.[Hl we have plotted Eq. dH) for i = 90° and i = 78°. For 
lower inclinations no eclipses occur for Ksac = 339.7 km s^^ 
and Porb = 0.134 d. For a given choice of the inclina- 
tion angle, Eq. ^ defines a unique relation Msec (Mwd)- 
If we further assume a mass for the secondary, we can 
investigate the possible range of the white dwarf mass 
for SDSS 0303+0054. For an extremely conservative lower 
limit of Msec ^ 0.08 Mp ( the lo wer mass limit for an 
M-dwarf, e.g. iDorman et al.l l|l989l) . with the spectrum of 
SDSS 0303+0054 clearly identifying the companion as a 
main-sequence star), Mwd ^ 0.68 M©. Even under this ex- 
treme assumption for the companion star, the white dwarf 
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Figure 8. Mass function plot (see Eq. |9j for SDSS 0303+0054, 
for i = 78° and i = 90°. Tlie gray liorizontal line corresponds 
to a lower mass limit for M-dwarfs. Tlie dotted horizontal line 
is the value of Msec, corresponding to a secondary spectral type 
of Sp(2) = M4, assuming a Sp(2) — Msec relation. The dashed 
vertical lines indicate the corresponding range for the white dwarf 
mass. 



in SDSS 0303+0054 has to be more massive than the av- 
erage field white dwarf (jKoester et al ] 1 19791 : iLiebert et al.l 
I2OO5I ). If we assume Sp(2) = M4, the upper limit of the 
spectral type according to our spectral decomposition of the 
SPSS spectrum, and use t he sp ectral type-mass relation of 
iRebassa-Mansergas et~all (|2007l '). we find Msec ^ 0.32 Mq, 
and hence Mwd ^ 0.96 Mq. These estimates assume i — 
90°, for lower inclinations the respective Mwd values be- 
come larger, as illustrated in Fig.[8l In summary, based on 
Kaec, Poib, and a generous range in possible A/scc we expect 
the white dwarf in SDSS 0303+0054 to be fairly massive, 
0.7Mq < Mwd < I.OMq. 



8 LIGHT CURVE MODEL FITTING 

Analysing the light curves of eclipsing binaries can provide 
strong constraints on the physical parameters of both stars. 
Here, we make use of a newly developed code written by 
one of us (TRM) for the general case of binaries contain- 
ing a white dwarf. The code offers the option to include 
accretion components (disc, bright spot) for the analysis of 
cataclysmic variables - given the detached nature of our 
targets, those components were not included. The program 
sub-divides each star into small elements with a geometry 
fixed by its radius as measured along the direction of centres 
towards the other star. The code allows for the distortion of 
Roche geometry and for irradiation of the main-sequence 
star using the approximation aTget. = aTstc + Fin, where 
Tgee is the modified temperature and Tsoc the temperature 
of the unirradiated companion, a is the Stefan-Boltzmann 
constant and Fi^r is the irradiating flux, accounting for the 



angle of incidence and distance from the white dwarf. The 
white dwarf is treated as a point source in this calculation 
and no backwarming of the white dwarf is included. The lat- 
ter is invariably negligible, while the former is an unneces- 
sary refinement given the approximation inherent in treating 
the irradiation in this simple manner. 

The code computes a model based on input system pa- 
rameters supplied by the user. Starting from this parameter 
set, model light curves are then fitted to the data using 
Levenberg-Marquardt minimisation, where the user has full 
flexibility as to which parameters will be optimised by the 
fit, and which ones will be kept fixed at the initial value. 

The physical parameters, which define the models, are 
the radii scaled by the binary separation, Rwo/a and 
Rscc/a, the orbital inclination, i, the unirradiated stel- 
lar temperatures of the white dwarf and the secondary 
star Tofi.wD and Tch, sec respectively, the mass ratio q — 
Msec /Mwd, To the time of mid-eclipse of the white dwarf 
and d the distance. We account for the distance simply 
as a scaling factor that can be calculated very rapidly for 
any given model, and so it does not enter the Levenberg- 
Marquardt optimisation. All other parameters can be ad- 
justed, i.e. be allowed to vary during the fit, but typi- 
cally the light-curve of a given system does not contain 
enough information to constrain all of them simultaneously. 
For instance, for systems with negligible irradiation, fitting 
Tcff.wD, T'eff.sec and d simultaneously is degenerate since a 
change in distance can be exactly compensated by changes 
in the temperatures. 

Some of our exposures were of a significant length com- 
pared to the length of the white dwarf's ingress and egress 
and therefore we sub-divided the exposures during the model 
calculation, trapezoidally averaging to obtain the estimated 
flux. 

Our aim is to combine the information from the ifscc 
radial velocity amplitudes determined in Sect.O the infor- 
mation contained in the hght curve, and the white dwarf 
effective temperature determined from the spectral fit in 
Sect.|S] to establish a full set of stellar parameters for the 
four eclipsing WDMS binaries. The adopted method is de- 
scribed in the next two subsections, where SDSS 0110+1326 
requires a slightly broader approach, as a correction needs to 
be applied to the Kscc velocity determined from the emission 
lines. 

8.1 SDSS 0303+0054, SDSS 1435+3733 and 
SDSS 1548+4057 

Our approach is to fit a hght curve model to the data for a 
selected grid of points in the Mwd — Msec plane. Each point 
in the A/wd — Msec plane defines a mass ratio q fEa. llO|l . and, 
using ifsec and Porb, a binary inclination i (Ea. lll|l . Hence, 
a light curve fit for a given (Mwd, Msec) combination will 
have Mwd, Msec, q, and i fixed, whereas Rwd, Race, Toff, sec, 
Tcff.wD, and To are in principle free parameters. In practice, 
we fix Teff^wD to the value determined from the spectro- 
scopic fit, which is sufficiently accurate, hence only Rwd, 
Rsac, Teff.sec, aud Tq are free parameters in the light curve 
fits. We leave To free to vary to account for the O — C errors 
in each individual light curve. The fitted values of Tq were 
of the order of 10 s, consistent with the O — C values quoted 
in Table[3] Each light curve fit in the Mwd — Msec plane 
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requires some initial estimates of -Rwd, Rsec, T^s, see- For 
-RwD we adopt the the oretical white dw arf mass-radius re- 
lation interpolated from lBergeron et a l.''s (199 5) tab l es. Fo r 
Reac we use the mass-radius relation of Baraf f e et alj (|l998l ). 
adopting an age of 5Gyr. For reff,scc, we use the spectral 
type of the secondary star as determined from the spectral 
decomposition (Sect.|6]l c ombined with the spectral type - 
temperature relation from Reba ssa-Mansergas et all (|2003)- 
-Rwd and Rscc are then scaled by the binary separation, 
Eq. 0. 

We defined large and densely covered grids of points in 
the MwD — Msec plane which generously bracket the initial 
estimates for Mwd and Msec based on the spectral decompo- 
sition/fitting. For SDSS 0303-1-0054, where no spectral fit to 
the white dwarf is available, we bracket the range in Mwd 
illustrated in Fig. [8] Points for which (formally) sini > 1 
were discarded from the grid, for all other points a light 
curve fit was carried out, recording the resulting Rwd, Race, 
and Teff.scc- The number of light curve fits performed for 
each system was between 7000 and 10000, depending on the 
system. 



8.2 SDSS 0110+1326 

In the case of SDSS 0110-1-1326 the inclination angle cannot 
readily be calculated through Eq. (|ll|l , as we measured K^ac 
from either the Ha or Ca ii emission lines, which do not trace 
the motion of the secondary's centre of mass, but the centre 
of light in the given emission line. The general approach in 
this case is to apply a correction to the observed Kscc, where 
the motion of the second ary star's centr e of m ass, K^ac .cor 
is expressed according to I Wade fc Horrid l| 19881 ) as 



Ks, 



Ks. 



1 - (1 + g) (Ai?/a) 



(13) 



where Ksec is the measured radial velocity, q is the mass 
ratio of the system, a the binary separation and AR is the 
displacement of the centre of light from the centre of mass of 
the secondary. AR can have a minimum value of zero, i.e. the 
two centres coincide and no correction is needed and a max- 
imum value of J?scc, i.e. all light comes from a small region 
of the secondary star closest to the primary. An assumption 
often used in the literature is that the emission due to irra- 
diation is distributed uniformly over the inner hemisphere 
of the secondary star, and is zero on its unirradiated side , 
in which case AR = (4/37r)flsec (e.g IWade fc Horndligsj : 
IWood et al.lll995bl : lOrosz et al.lll999l : IVennes et al.lll999l ). A 
more physical model of the distribution of the irradiation- 
induced emission line flux can be derived from the analysis of 
the orbital variation of the equivalent width of the emission 
line. A fine example of this approach is the detailed study of 
the Ca 11 AA 8498.02, 8542.09,86 62.14 emissiori i n the WDMS 
binary HR Cam presented by iMaxted et al.l |l99^) . How- 
ever, given the small number of spectra, it is currently not 
possible to apply this method to SDSS 0110-1-1326. 

Our approach for modelling SDSS 0110-1-1326 was to 
calculate the binary separation from Eq. (|SJ for each pair of 
(Mwd, Mace) of the initial grid and then, using Eq. (|13p . cal- 
culate various possible -K'sec.cor values for this (Mwd, Msec) 
point, assuming various forms of AR. We adopted the follow- 
ing three different cases: (i) AR = 0, i.e. the centre of light 
coincides with the centre of mass, (ii) A_R = (4/37r)_Rsec, 
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Figure 9. Grid of (A/wDiMsec) points for which the light curve 
of SDSS 0110+1326 has been fitted. Different assumptions for the 
i^-correction were made: AR = (black strip), AR = (4/37r)/?soc 
(dark gray strip) and AR = Rsec (light gray strip). Loft panel: 
(Mwd, Msec) grid for Kscc = KsecHa = 200.1kms~l. Right 
panel: (Mwdi^-^scc) grid for Kscc = ^^scc.Call = 178.8kms~l. 



Table 6. Corrected values of the radial velocity of the secondary 
star in SDSS 0110+1326, after assuming a form of AR, for both 
the cases of the Ha and the Can emission lines. 



AR 


mean [kms ^] 


std.dev. 


range [kms ^] 


Ha 











200.1 






(4/37r) Rsec 


213.1 


2.7 


207.9-218.2 


Rscc 


225.2 


3.5 


218.5-232.1 


Can 











178.8 






(4/37r) Rsec 


195.1 


4.2 


186.6-202.7 


Rsec 


209.1 


6.1 


197.2-219.9 



i.e. the uniform distribution case and (iii) AR =i?sec, i.e. 
the maximum possible displacement of the centre of light. 
This was done for both the Ha and the Can lines. Having 
attributed a set of corrected radial velocity values to each 
point of the grid, we could now make use of Eq. (Illf) to cal- 
culate the corresponding inclination angles. Again, points 
with a (formal) sin i > 1 were eliminated from the grid. The 
allowed (MwD,Afscc) pairs for both lines and all three AR 
cases are shown in Fig.|9] 

The first of the three AR cases, AR — 0, is obviously 
not a physically correct approach, as it implies that no cor- 
rection is necessary. However, we did use it as a strict lower 
limit of the radial velocity of the secondary. In this sense, the 
third case, AR =Rsec, is a strict upper limit for Ksec- Ta- 
blelUlists the ranges Ksec,coT corresponding the our adopted 
grid in the (Mwd, Afsce) plane shown in Fig. [9] 
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9 RESULTS 

Our light curve fitting procedure outfined in Sect .[8] yields 
fitted values for Rwu, -Rscc, and Toff, sec for a large grid in the 
M^jYY) — Mace plane, where each {Mwu, Msec) point defines 
q and i. Analysing these results was done in the following 
fashion. 

In a first step, we applied a cut in the quality of the 
light curve fits, where we considered three different degrees 
of "strictness" by culling all models whose was more than 
1, 2 or 3(T above the best-fit value. This typically left us with 
a relatively large degeneracy in Mwd and Msec . 

In a second step, we had to select among various equally 
good light curve fits those that are physically plausible. 
Given that our set of fitted parameters was underconstrained 
by the number of observables, we had to seek input from the- 
ory in the form of m ass-radius relations fo r the white dwarf 
(using the tables of iBergeron et al. || 19951) and for th e M- 
dwarf (using the 5.0 Gvr model of lBaraffe et all 119981 ). We 
then calculated for each model the relative difference be- 
tween the theoretical value of the radius and the value of 
the radius from the fit, 



0.5 



5R = 



Rat — Rth 



Rth 



(14) 



where Rfit and Rth are the radius obtained from the light 
curve fit, and the radius from the mass-radius relation, re- 
spectively. Assuming that the binary components obey, at 
least to some extent, a theoretical mass-radius relation, we 
defined cut-off values of SR = 5%, 10%, and 15%, above 
which light curve models will be culled. The maximum of 
SR = 15% was motivated by the radius excess over theoret- 
ical main-sequence models observed in eclipsing low-mass 
binaries (Ribas ct al. 2007), and by the current observa- 
tional constraints on t he white dwarf mass-radius relation 
jProvencal et al.lll998l ). This cut in radius was applied in- 
dividually for the white dwarf, (SRwd) and the companion 
star (SRscc)- 

Combining both the ^-nd SR cuts left us in gen- 
eral with a relatively narrow range of system parameters 
that simultaneously satisfy the radial velocity amplitude 
and the morphology of the light curve. In a final step, we 
examined how well the stellar masses determined from the 
light curve/radial velocity analysis agreed with those derived 
from the spectral decomposition. 



9.1 SDSS 0303+0054 

Modelling the light curve of SDSS 0303-1-0054 involved the 
following two problems. Firstly, the time resolution of our 
data set poorly resolves the white dwarf ingress and egress 
phases, and consequently the white dwarf radius can only be 
loosely constrained. We decided therefore to fix the white 
dwarf radius in the light curve fits to the value calcu- 
lated from the adopt ed white dwarf mass-radius relation of 
IBergeron et al.l Il995l . Secondly, the light curve displays a 
strong, but slightly asymmetric ellipsoidal modulation, with 
the two maxima being of unequal brightness (Fig.[l]). Similar 
light curve morpholo gies have been observed in the WDMS 
binaries BPM 71214 ^Kawka fc Vennesll2003D and LTT560 
jTappert et al.|[2007l '). and have been attributed to star spots 
on the secondary star. By design, our light curve model can 
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Figure 10. Liglit curve model fitting results for SDSS 0303-f 0054. 
Left panel: Mwd and Mgec values corresponding to fits with 
values within Itr of the minimum value (gray points) and, 
simultaneously, with SRsec ^ 0.05 (black points). Right panel: the 
same, only in the q — i plane. Also depicted in the left panel are 
curves corresponding to the mass function (solid black lines, i = 
90° and 78°) which (by definition) bracket the possible solutions, 
Sp(2) — Msec relations (dotted, horizontal, black lines) and the 
range of possible (AfwDiA/scc) values (dashed, horizontal and 
vertical, black lines) 



not provide a fit that reproduces the observed asymmetry. 
Having said this, the presence of ellipsoidal modulation pro- 
vides an additional constraint on R^cc that is exploited by 
the light curve fit. 

A number of models passed the strictest configuration 
of our cut-offs, i.e. la and SR^cc = 5% (no cut was used in 
SRwD, as i?wD has been fixed during the fits). The possible 
range of solutions in the Mwd — Msec plane is shown in 
Fig.M 

Gray dots designate those light curve fits making the 
la cut, black dots those that satisfy both the la and 
SRscc ~ 5% cuts. The resulting ranges in white dwarf masses 
and secondary star masses are MwD.icflt = 0.88 — 0.95 M© 
and Msec = 0.22 — 0.28 M©, respectively, corresponding to 
a white dwarf radius of Rwu = 0.008 — 0.009 R© and a 
secondary radius of Rscc = 0.25 — 0.27 R©. Fig. 1111 shows 
one example of the light curve fits within this range for 
the model parameters A4"wd = 0.91 M©, i?wD ~ 0.009 Rq, 
Msec = 0.25 Mq, Rsec = 0.26 R© and i = 82.6°. 

To check the consistency of our light curve fits with 
the results from the spectral decomposition, we indicate in 
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Figure 11. Model fit to the CA2.2 filterless light curve of 
SDSS 0303+0054, for A/wd = 0.91 Mq, Rwd = 0.009 Rq, 
Msec = 0.25 M0, i?scc = 0.26 R© and i = 82.6°. The model 
meets both our (within la) and the SR (within 5%) cut-offs. 
The residuals from the fit are shown at the bottom of the panel. 
Inset panel: data points and model fit focused around the eclipse. 



Fig. [To] the radii of M-dwarfs with spectral types Sp(2) = 
M4 — M6 in steps of 0.5, based on the spectral type- 
mass relation given bv lRebassa-Mansergas et al.l l|2007h . The 
secondary star mass from the hght curve fit, Msec.icfit = 
0.23 — 0.28 M0, corresponds to a secondary spectral type 
of Sp(2)icflt = M4.5 — M5.5. The expected spectral type of 
the secondary from the spectral decomposition (Sect.O was 
SPapflt(2) = M4.5±0.5, hence, both methods yield consistent 
results. 



9.2 SDSS 1435+3733 

For SDSS 1435+3733, which is partially eclipsing, our tem- 
poral resolution, although again not ideal, was nevertheless 
deemed to be adequate to constrain the white dwarf radius 
from the light curve fits. Thus, a SRwu cut was also applied. 

We applied again our strictest cuts on the light curve 
models, la and SRwd = SRsec = 5%, and the parameters of 
the surviving models are shown in Fig. 1121 where the mean- 
ing of the symbols is the same as in Fig. 1101 

These fits imply a white dwarf mass and radius of 
MwD = 0.48-0.53 Mq and i?wD = 0.014-0.015 R©, respec- 
tively, and a secondary star mass and radius of Msec = 0.19— 
0.25 M© and Rscc~ 0.22 — 0.25 R©, respectively. A sample 
fit that obeyed all three constraints is shown in Fig. 1131 The 
model parameters are Mwd ~ 0.5 M©, Rwu = 0.015 R©, 
Msec = 0.21 Mq, i?scc= 0.23 R© and i = 77.6°. 

Comparing the white dwarf masses from the spectro- 
scopic decomposition/fit, MwD.spflt = 0.41 ± 0.08 Mq, with 
the range of white dwarf masses allowed by the light curve 
fitting, MwD.icflt ~ 0.48 — 0.53 Mq, reveals a reasonable 
agreement. Regarding the spectral type of the secondary 
star, we indicate in Fig. [12] again the masses of M-dwarfs in 
the range M4-M6, in steps of 0. 5 spectral classes, follow- 
ing th e Sp(2) — Msec relation of iRebassa-Mansergas et ah] 
(|2007l ). This illustrates that the light curve fits result in a 
Sp(2)icfit ~ M5 — M6, whereas the spectral decomposition 
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Figure 12. Light curve model fitting results for SDSS 1435+3733. 
Left panel: Mwd and Msec values corresponding to fits with 
values within Itr of the minimum value (gray points) and, 
simultaneously for both white dwarf and secondary radii, with 
SR ^ 0.05 (black points). Right panel: the same, only in the q — i 
plane. Also depicted in the left panel are mass functions (i = 70° 
and i = 80°, solid black lines), Sp(2) — Mace relations (dotted, 
horizontal, black lines) and the range of possible (Mwdi-'^soc) 
values (dashed, horizontal and vertical, black lines) 



provided Sp(2)spfit = M4.5 ± 0.5. Hence, also the results for 
the secondary star appear broadly consistent. 

ISteinfadt et al.] (|2008h first reported the eclipsing na- 
ture of SDSS 1435+3733, and estimated Mwd = 0.35 - 
0.58 Mq, J?wd = 0.0132-0.0178 Rq, Afscc = 0.15-0.35 Mq 
and _Rscc = 0.17—0.32 Rq. The parameter ranges determined 
from our analysis (Mwd = 0.48 - 0.53 Mq, Rwt> = 0.014 - 
0.015 Rq, Msec = 0.19- 0.25 Mq and fise c= 0.22- 0.25 Rq) 



are fully consistent with ISteinfadt et al.'s work. In Fi g. [T4l 
we overplot on the photometrv of Isteinfadt et al.l (l2008l ) the 
light curve model shown in Fig. [13] along with our data, il- 
lustrating that our solution is consistent with their data. 

While the te mporal resolution of our photometry is 
worse than that o f Isteinfadt erall (|2008l ). our analysis bene- 
fitted from two additional constraints, firstly the mass func- 
tion (Eq.[9]) determined from our spectroscopy, and sec- 
ondly the detection of a weak ellipsoidal modulation in 
ou r J-band light c urve (Fig.ll3 p . While the observations 
of Isteinfadt et al.l l|2008l ) covered the entire binary orbit, 
their data were obtained through a BG-39 filter, centred 
at 4800 A, where the fiux contribution of the companion 
star is negligible. A well-established difficulty in modelling 
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Figure 13. Model fit to IAC80 /-band light curve of 
SDSS 1435+3733, for Mwd = 0.5 M©, RwD = 0.015 R©, Msec = 
0.21 M0, Raec = 0.23 R0 and i = 77.6°. The model meets both 
the Icr cut in x'^, and a SR sC 0.05 cut for both the white dwarf 
and the secondary star. The residuals from the fit are shown at 
the bottom of the panel. Inset panel: data points and model fit 
focused around the eclipse. 
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Figure 15. Light curve model fitting results for SDSS 1548+4057. 
Left panel: Mwd and Mgec values corresponding to fits with 
values within Itr of the minimum value (gray points) and, 
simultaneously for both white dwarf and secondary radii, with 
SR ^0.1 (black points). Right panel: the same, only in the q — i 
plane. Also depicted in the left panel are mass function for i = 78° 
and i = 90° (solid black linos), Sp(2) — Msec relations for spectral 
types M6, 6.5, and 7 (dotted, horizontal, black lines) and the 
range of possible (Mwdi-'^soc) values (dashed, horizontal and 
vertical, black lines). 
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Figure 14. Our light curve model for SDSS 1435+3733, as de- 
termined from fitting our IAC80 J-band photometry a s shown in 
Fig-HSl along with the data from lSteinfadt et al]||2008^ . illustrat- 
ing a high degree of consistency. 

light curves of partially eclipsing binary stars is the fact 
that, although the sum of the two radii can be accurately 
defined, their ratio ca n only be loosely constrained (e.g. 
ISouthworth et alj|2007h . 7?WD and Rscc are strongly corre- 
lated, and fitting both of them leads to degeneracy. Addi- 
tional information is needed to lift this degeneracy, which 
is mainly pr ovided by ellipsoidal modulation or reflection 
effects (e.g. iHilditch et al..,1996a ibil. 

9.3 SDSS 1548+4057 

For SDSS 1548+4057, we follow a very similar approach as 
for SDSS 1435+3733, adopting a la cut in but a slightly 
less strict cut on the radii, SRwd ~ SRscc ~ 10%. FigurefTSl 
shows the solutions of our hght curve fits that survived those 



criteria, where plot symbols have the same meaning as in 
FigsfT0land[T2] 

The models satisfying both cuts imply a white dwarf 
mass and radius of Mwd = 0.61 — 0.68 Mq and Rwd ~ 
0.010 — 0.011 Rq, respectively, and a secondary mass and 
radius of M^ec = 0.15- 0.20 Mq and 7£sec = 0.17- 0.20 Rq, 
respectively. A sample fit, that obeyed all three constraints, 
is shown in Fig. [TB] The model parameters are Mwd ~ 

0.65 Mq, RwD = 0.011 Rq, Msec = 0.18 Mq, i?sec = 0.18 Rq 

and i = 85.6° . 

The white dwarf parameters obtained from the light 
curve fitting are in full agreement with those derived from 
the spectral decomposition. The former give a possible range 
of white dwarf masses of MwD.icflt = 0.61 — 0.68 Mq, while 
the latter yield a value of MwD.spfit = 0.62 ± 0.28 Mq. 
In Fig. 1151 we indicate the masses corresponding to sec- 
ondary spectral types M6, M6.5, and M7. The spectral fit- 
ting points to a Sp(2)apflt ~ M6 secondary. Our models in- 
dicate Sp(2)icflt = M6 — M7 again consistent with the spec- 
troscopic result. 
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Figure 16. Model fit to the WHT/iJ-band fight curve of 
SDSS 1548+4057, for A/wd = 0.65 Mq, i?wD = 0.011 R©, 
Msec = O.ISMq, Rscc = O.183R0 and i = 86.6°. The model 
meets both a Icr cut-off, and the &R 0.10, cut-off, for both 
the white dwarf and the secondary radii. The residuals from the 
fit are shown at the bottom of the panel. Inset panel: data points 
and model fit focused around the eclipse. 



Table 7. Summary of the fitting resuhs of SDSS 0110-f 1326, for 
all three AR prescriptions, for both the cases of the H« and the 
Call emission lines. In all cases we applied Scr and &R = 15% 
cut-offs. A {^J) means some models passed the cut-off, whereas a 
(-) means no models passed. 
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Figure 17. Results of the model fitting of SDSS 0110-f 1326 for 
Ha (left panels) and Call (right panels), for a Scr confidence level 
and 5R < 0.15. Black points correspond to AR = 0, gray to 
AR = (4/37r)i?scc (where applicable, see text for details). Top 
panels show the (M^Di Msec) plain, bottom panels the q — i plain. 



9.4 SDSS 0110+1326 

Because of the uncertainty in the true Kaec velocity of the 
secondary star, we inspected a total of six different grids of 
light curve fits for SDSS 0110-1-1326 (three K-corrections for 
each of the observed Kscc, Ha and Kaec, Caii, see Sect. l8.2l and 
TableO, and consequently, this star needs a slightly more 
extensive discussion compared to the other three systems. 

In a first exploratory step, we applied our least strict 
cut-offs, adopting a 3a cut on the fit quality, and a 5R = 
15% cut for both the white dwarf and the secondary star 
radii. Even with these loose constraints, the AR = i?scc 
case, i.e. assuming that the emission is concentrated on the 
surface of the secondary closest to the white dwarf, does not 
produce any solutions that meet the cut-off criteria. For the 
uniform illumination case, AR — (4/37r) f?scc, there were no 
solutions for J^soc ~ -K^socHa, but a few solutions for the 
T^sec = KsBc,Ca,u- Finally, a number of solutions existed for 
the case AR = 0, i.e. assuming no K-correction, for both Ha 
and Can. Table[7] summarises this analysis, and the models 
that survived the culling defined ranges of (MwD,Afscc) pairs 
which are illustrated in Fig. 1171 

Given that the uniform irradiation prescription, AR = 
(4/37r) J?sec represents a mid-way between the two extreme 
assumptions AR — i?scc and AR — 0, and that a reason- 



ably large range of possible solutions was found for the Ca ii 
velocities (Fig. I17|l . we decided to explore tighter con- 
straints for this case. While a number of light curve models 
are found within la of the minimum value, both SRwd 
and SRbcc were ~ 15% in all cases. No models survived if 
we reduced the value of SR in addition to the tighter 
constraint. 

As an illustration of the quality of the light curve fits 
achieved. Fig. 1181 shows a fit with Mwd = 0.59 Mq, Rwd = 
0.015 R©, Msec = 0.28 M©, TJsec = 0.32 R©, i = 85° and 
A'scc = 197 km s-^ The mass of the secondary corresponds 
to a spectral type of Sp(2) — M4.5. We stress, however, 
that we do not adopt this model as a final solution for the 
parameters of SDSS 01104-1326, but rather consider it as a 
physically plausible set of parameters. 

Taking the la confidence levels of the light curve fits 
at face value, they imply a white dwarf mass of MwD.icflt = 
0.59 — 0.7 M© and a secondary mass of Mscc,icflt ~ 0.28 — 
0.38M©. 

For the white dwarf, these numbers are in contrast with 
the results from the spectral decomposition, MwD,spflt = 
0.47 ± 0.02 M©. The errors in the spectroscopic white dwarf 
parameters are purely of statistical nature. However, even 
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Figure 18. Model fit to the CA2.2 /-band light curve of 
SDSS 0110+1326, for A/wd = 0.59 Mg, Rwd = 0.015 Rq, 
Msec = 0.28 M0, iJsoc = 0.32 R0 i = 85°, and i^sec = 
197kms~^. The model radii agree with the theoretically pre- 
dicted ones within 15%, for both the white dwarf and the sec- 
ondary star. The residuals from the fit are shown at the bottom 
of the panel. Inset panel: data points and model fit focused around 
the eclipse. 



assuming a reasonably large systematic error of 0.2 dex in 
loggr (implying MwD.spfit = 0.39 — 0.55 Mq) the white dwarf 
masses implied by the light curve fit and the spectroscopic 
fit appear to be inconsistent. 

Regarding the secondary, adopting A/scc,icflt = 
0.28 - 0.38 Mq and the S p(2) - Msec relation of 
iRebassa-Mansergas et al]|2007l suggests a spectral type of 
Sp(2)icflt = M3 — M4.5, which is consistent with the results 
of the spectral decomposition, Sp(2)spfit — M4 ± 1. 

While our analysis of the light curve fits demonstrates 
that the data can be modelled well, it is clear that the lack of 
a reliable value for Kscc prevents the choice of the physically 
most meaningful solution. Consequently, a definite determi- 
nation of the stellar parameters of SDSS 0110-1-1326 is not 
possible on the basis of the currently available data, and we 
adopt as a first cautious estimate the values of our spectral 
decomposition fit. As mentioned in Sect. 18.21 a physically 
motivated A'-correction can be modelled once a larger set of 
spectra covering the entire binary orbit are available. 



10 POST COMMON ENVELOPE EVOLUTION 

Considering their short orbital periods, the four new eclips- 
ing WDMS binaries discussed in this paper must have 
form e d through comm on envelope evolution l|Paczvnskil 
1 19761 : IWebbinkI |20o3). Post common envelope binaries 
(PCEBs) consisting of a white dwarf and a main sequence 
star evolve towards shorter orbital periods by angular mo- 
mentum loss through gravitational radiation and magnetic 
braking until they enter the semi- detached cataclysmic 
variable configuration. As shown by ISchreiber fc Gansickd 
l|2003t ). if the binary and stellar parameters are known, it is 
possible to reconstruct the past and predict the future evolu- 
tion of PCEBs for a given angular momentum loss prescrip- 
tion. Here, we assume c lassical disrupted magnetic braking 
l|Rappaport et al]ll983l ). i.e. magnetic braking is supposed 



Table 8. Reconstructed and predicted orbital periods as well as 
expected evolutionary times for the four new eclipsing WDMS 
binaries. The cooling age of the white dwarfs is icool while Pce 
and Pgtj are the orbital periods at the end of the common envelope 
phase and at the onset of mass transfer, respectively. The time 
still needed to enter the semi-detached CV configuration is tgj- 



SDSS J 


^cool [yr] 


Psd [d] 


Pce [d] 


tsd [yr] 


0110-1-1326 


1.41e+07 


0.097 


0.332 


1.78e-|-10 


1548-1-4057 


3.91e+08 


0.069 


0.191 


4.23e-|-09 


0303-1-0050 


3.23e+09 


0.097 


0.228 


5.95e-|-08 


1435-f3733 


2.53C+08 


0.088 


0.132 


9.50e-|-08 



to be much more efficient than gravitational radiation but 
only present as long as the secondary star has a radiative 
core. For fully convective stars gravitational radiation is as- 
sumed to be the only angular momentum loss mechanism, 
which applies to all four systems. The current ages of the 
PCE Bs are determined b y interpolating the cooling tracks 
from I Wood et all (|l995a'). Table|8] lists the evolutionary pa- 
rameters our four eclipsing WDMS binaries. 

The calculated evolutionary tracks of the four systems 
are shown in Fig. 1191 The present position of each binary is 
indicated by an asterisk. The grey-shaded region indicates 
the 2 — 3 h orbital period gap, i.e. the orbital period range 
where only a small number of cataclysmic variables (CVs) 
has been found. Given that the four WDMS binaries contain 
low-mass secondary stars, all of them are expected to start 
mass transfer in or belo w the period gap. Thi s unde rlines 
the fact first noticed by ISchreiber fc Gansickd (J2003) that 
only very few progenitors of long-period CVs are known. 
SDSS 0303-1-0050 has passed most of its PCEB lifetime and 
its current orbital period differs significantly from the recon- 
structed value at the end of the CE phase. The other three 
systems are rather young PCEBs with a current orbital pe- 
riod close to Pce- 



11 CONCLUSIONS 



Since iNelson fc Youn j l|l970h reported the discovery of 
V471 Tan, as the first eclipsing white dwarf main se- 
quence binary, only 7 additional similar systems have 
bee n discovered, t h e late st being SDSS 1435-1-3733 found 
by ISteinfadt et al.l (|2008l ) . Here, we presented follow-up 
observations of SDSS 1435-1-3733, which we independently 
identified as an eclipsing WDMS binary, and of three 
new discoveries: SDSS 0110+1326, SDSS 0303+0050, and 
SDSS 154 8+4057. A fifth system has just be en announced by 
our team (|Nebot Gomez-Moran et~al ] |2008l). and two eclips - 
ing WDMS binaries were identified bv iDrake et all (|2008l ) 
in the Catalina Real-Time Transient Survey. Table[5] lists 
the basic physical parameters of the fourteen eclipsing white 
dwarf main sequence binaries currently known. 

Inspection of Table|9] shows that the orbital periods of 
the known eclipsing binaries range from three to 12 hours, 
with three of the systems presented here settling in at the 
short-period end. 

The white dwarf masses cover a range ~ 0.44 — 0.9 Mq, 
with an average of < AfwD >= 0.57±0.16 Mq, which is only 
slightly lower than the average mass of single white dwarfs 
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Table 9. Physical parameters of the known ecHpsing white dwarf main sequence binaries. A dash sign means that no corresponding value was 
given by the respective authors. Errors on the parameters are also shown, if quoted by the authors. 
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CSS080408t 




0.40 + 0.05 


0.022 + 0.003 


32595 + 675 


7.36 + 0.17 


0.26 


0.27 





Kacc [kms 



-11 



Sp2 



150.7+1.2 
148.46+0.56 
136.4+0.8 
166+4 



181 + 3 
180-200 

195.8* 

190+9 
190.2+3.5 

274.7 
241+8.1 

266+5 

339.7 
310+10 
289.3+12.9 



260.9 



K2 
K2 
M3.5-M4 

M3V 
M0-M6 
M3-M5 
M4+1 
M3-M5 
M6 
M6 
M4 
M5.5-M6.5 
M2-M4 
M3.5-M4 

M4 
M4-M5 
M4.7-M6.1 
M4.5-M5 
M4-M6 
M4-M5 
M5 



Ref erences: (lallYoung fc NeIsonll972l: f Ib'llO'Brien et al.l200ll; (2') lMaxted et al.l2004l: (3)\van den Besselaar et al.l2007l; (4a ) lGreen et al.ll97L. 
f4b')|FuIbright et al.lll993l: fSllNebot Gomez-Moran et a.lJ l2008l: ( 6) this work; (7a) Bruch fc Diaj Il998l: CTbteruchl Il999l: f7c')lMaxted et aL 
I2OO7I : (Safcawka et al.l l2002ir78b')IO'Donoghue et al.l |2003|: (gl lOrake et al.l bOOsI: (lOal lWood fc MarshI Il99ll : (lObl lCatalan et al.l Il994l : 
flOc^ lHaefner et al .I I2OO4I: flD lSteinfadt et al.ll2008l 

Notes: * Calculated from Ksini provided by the authors; + CSS080502:090812+060421 = SDSS J090812.04+060421.2, we de- 
termined the white dwarf and secondary star paremeters from decomposing/fitting the SDSS spectrum (see Sect.|6] and 
iRebassa-Mansereas et al] l2007t ). these parameters should be taken a s guidance as no light c u rve m odelling has been done so far; t 
CSS080408: 142355+240925 = SDSS J142355.06+240924.3, also found bv iRebassa-Mansergas et al.l ll2008al '). white dwarf and secondary star 
parameters determined as for CSS080502. 




Figure 19. Evolutionary post common envelope tracks for our 
four eclipsing binaries. From left to right: SDSS 1435+3733, 
SDSS 1548+4057, SDSS 0303+0054 and SDSS 0110+1326. The 
current position of the systems is indicated by an asterisk. We 
assumed only gravitational wave radiation as angular momentum 
loss agent as all four systems are thought to contain fully convec- 
tive secondary stars. The orbital period gap of CVs is highlighted 
in gray - all four systems will start mass transfer in, or shortly 
below the gap. 



(|Liebert et al.ll2005l : iFinlev et al]| 19971 : iKoester et aLlll979l ) . 

This finding is somewhat surprising, as all these eclipsing 
WDMS binaries have undergone a common envelope evo- 
lution, which is expected to cut short the evolution of the 
primary star and thereby to produce a significant number of 
low-mass He-core white dwarfs. Our estimates for the white 
dwarf masses in SDSS 0110+1326 and SDSS 1435+3733 sug- 
gest that they contain He-core white dwarfs. The white 
dwarf in SDSS 1548+4057 is consistent with the average 
mass of single white dwarfs. SDSS 0303+0054 contains a 
fairly massive white dwarf, potentially the most massive one 
in an eclipsing WDMS binary. 

The secondary star masses are concentrated at very 
low masses, the majority of the fourteen systems in Ta- 
ble|§]have A'Isac < 0.6 M©, and ~ 9 systems (including the 
four stars discussed here) have Msec ^ 0.3 Mq. This is a 
range w here very few low-mass stars in eclipsing binaries are 
known (lRibas|[2006l ). making the secondary stars in eclips- 
ing WDMS binaries very valuable candidates for filling in 
the empirical mass-radius relation of low mass stars. This 
aspect will be discussed in a forthcoming paper. 

The distribution of secondary star masses among the 
eclipsing PCEBs in Table [9] is subject to similar se- 
lection effects as th e who le sample of known PCEBs. 
ISchreiber fc Gansickd l|2003D found a strong bias towards 
late-type secondary stars in the pre-SDSS PCEB pop- 
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ulation, which arose as most of those PCEBs were 
identified in blue-colour surveys. T he currently e merg- 
ing population of SP SS PCEBs ()Silvestri et all l2006l : 
iRebassa-Mansergas et al. 2008a ) should be less biased as the 
ugriz - colour-space allows to identify WDMS binaries with 
a wide range of secondary spectral types and white dwarf 
temperatures. 

Future work will need to be carried out on two fronts. 
On one hand, more accurate determination of the masses 
and radii will have to be made to unlock the potential 
that WDMS binaries hold for testing/constraining the mass- 
radius relation of both white dwarfs and low mass stars. The 
key to this improvement are high-quality light curves that 
fully resolve the white dwarf ingress/egress, and, if possible, 
the secondary eclipse of the M-dwarf by the white dwarf. 
In addition, accurate Kwo velocities are needed, which will 
be most reliably obtained from ultraviolet intermediate res- 
olution spectroscopy. On the other hand, it is clear that the 
potential of SDSS in leading to the identification of addi- 
tional eclipsing WDMS binaries has not been exhausted. Ad- 
ditional systems will be identified by the forthcoming time- 
domain surveys, and hence growing the class of eclipsing 
WDMS binaries to several dozen seems entirely feasible. 
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